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We study the black hole particle production in a regular spacetime metric obtained in a minisu- 
perspace approach to loop quantum gravity. In different previous papers the static solution was 
obtained and shown to be singularity-free and self-dual. In this paper expanding a previous study 
of the black hole dynamics we repeat the Hawking analysis which leads to a thermal flux of particles 
at the future infinity. The evaporation time is infinite and the unitarity is recovered due to the 
regularity of the spacetime and to the characteristic behavior of the surface gravity. 



I. INTRODUCTION 

Black holes are one of the most fascinating predictions 
of Einsten's gravitational theory. Today we know that 
they are not just a mathematically possible solution of 
general relativity, but part of Nature. Since more than a 
decade now, we have good evidence that our Milky Way, 
as other galaxies, hosts many stellar black holes as well 
as a supermassive black hole in its center. 

From the perspective of quantum gravity, black holes 
are of interest because of the infinite curvature towards 
their center which signals (probably) a breakdown of 
General Relativity. It is an area where quantum grav- 
ity effects are strong, and it is generally expected that 
these prevent the formation of the singularity. 

In the mid-seventies Steven Hawking [l[ and Jacob 
Bekenstein Q showed that since classical black holes ra- 
diate particles, their mass slowly decrease and finally dis- 
appear, leaving a paradox. Supposing that matter which 
enters the black hole was in a pure quantum state at early 
time, then the existence of an horizon transforms such 
matter into the mix thermal state of Hawking radiation 
at late time for an external observer, and the complete 
evaporation of the singularity destroys the correlations 
between the radiation and the information "swallowed" 
by the black hole that would allow to reconstruct the 
original pure quantum state. The role of the singularity 
is crucial in this paradox. Suppose there is an entan- 
gled pure state at early time, and a part of the entangled 
system ends inside the black hole while the other part 
remains outside the event horizon, the result is a mixed 
state after the partial trace is taken over the interior of 
the black hole. Now, since everything that falls inside 
the black hole reaches the singularity (r = 0) in finite 
proper time, the part of the Hilbert space that is traced 
over disappears and never appears again with the disap- 
pearance of the black hole. In the standard semiclassical 
treatment this is the scenario, since the black hole emits 
particles in the form of Hawking radiation and the hori- 
zon radius decreases and approaches the singularity until 
both, singularity and horizon, vanish in the endpoint of 
evaporation Q. However, if the singularity does not ex- 



ist, this scenario cannot be correct. Since the singularity 
plays a central role for the causal spacetime diagram, its 
absence in the presence of quantum gravitational effects 
has consequences for the entire global structure and 
its removal is essential for resolving the black hole infor- 
mation loss problem 0, @ . To understand the dynamics 
of the gravitational and matter fields, it is then necessary 
to have a concrete model. 

It is thus promising that a resolution of the big bang 
as well as the black hole singularites [§4ic| has been 
achieved in a simplified version of loop quantum gravity 
LQG) [HI, known as loop quantum cosmology (LQC) 
The regular static black hole metric was recently 
derived in [H | 12 1 , and studied more closely in [l3[ . The 
polymeric approach to solve the black hole singularity 
problem was also applied in [l4j]. An alternative reso- 
lution of the black hole singularity was obtained in an 
effective, noncommutative approach to quantum gravity 
[Tol l and in asymptotically safe quantum gravity jl6| . In 
others works [6j [171 ] , a 2-dimensional model was used to 
study the evaporation process in the absence of a sin- 
gularity. Recently also regular spinning loop black holes 
were obtained in [lj| ■ 

Here, we will use a 4-dimensional model based on the 
static solution derived in [ll[ and generalize it to a dy- 
namical case which then allows us to examine the causal 
structure. This generalization holds to good accuracy in 
all realistic scenarios. This approach should be under- 
stood not as an exact solution to a problem that requires 
knowledge of a full theory of quantum gravity, but as a 
plausible model based on preliminary studies that allows 
us to investigate the general features of such regular black 
hole solutions with quantum gravity corrections inspired 
by LQG. The main goal of this paper is to reproduce 
the Hawking calculation of particle creation in a partic- 
ular regular black hole background and to show that the 
whole process, collapse and complete evaporation, is uni- 
tary. In this analysis we will make extensively use of the 
Fabbri and Navarro-Salas book [2(j ■ 

Non-singular black holes were considered already by 
Bardeen in the late 60s and have a long history [l9|, l2lT - 
|3(I In this paper we extend the procedure in 
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[301 ] [3l| analyzing the particle production and the uni- 
tary problem as result of the complete evaporation. The 
paper is organized as follows. We start in the section ITT1 
by recalling the regular static metric we will use in the 
whole article. In section IIIII we generalize it to a col- 
lapse scenario as already done in [3l| . In section IIVI we 
derive and solve the scalar field equations of motion in 
this regular black hole background and in section [V] wc 
use these solutions to reproduce the Hawking analysis 
leading to a thermal spectrum. In section IVII we sum- 
marize the complete dynamics for the LBH collapse and 
evaporation. We dedicate section [VTT] to review the basic 
aspects of the unitarity and information loss problems 
associated to classical black holes and finally in section 
I Villi we show how in this quantum corrected model these 
problems are solved. The signature of the metric used in 
this paper is ( — ,+,+,+) and we use the natural unit 
convention h = c = Gn = 1. 



the fixed graph decided above. The diff-constraint 
is then identically zero because of homogeneity and 
the Gauss constraint is zero for the Kantowski- 
Sachs spacetime. 

ii We solve the Hamilton equation of motion for 
the holonomic Hamiltonian system imposing the 
Hamiltonian constraint to be zero. 

iii The third step consists to extend the solution to 
the whole spacetime; this step is mathematically 
correct but, since we found the solution in the ho- 
mogeneous region, we can not be sure of the cor- 
rect Hamiltonian constraint polymerization in the 
full spacetime. However we believe such polymer- 
ization exists. 

This quantum gravitationally corrected Schwarzschild 
metric can be expressed in the form 



II. THE REGULAR BLACK HOLE SPACETIME 

LQG is a candidate theory of quantum gravity. It is 
obtained from the canonical quantization of the Einstein 
equations written in terms of the Ashtekar variables [37| , 
that is in terms of an su(2) 3-dimensional connection A 
and a triad E. The result [39| is that the basis states of 
LQG are closed graphs made of edges and nodes labelled 
by irreducible SU(2) representations and intcrtwiners re- 
spectively. Physic ally, the edges represent quanta of area 
with area 87T7Zp y/ j{j + 1), where j is the representation 
label of the edge (an half- integer) , lp is the Planck length, 
and 7 is a parameter of order 1 called the Immirzi param- 
eter. Vertices of the graph represent quanta of 3- volume. 
The area is quantized and the smallest possible quanta 
correspond to an area 



^min = 47ta/3Vp ■ 



(1) 



To obtain the simplified black hole model, the following 
assumptions were made. First, the number of variables 
was reduced by assuming spherical symmetry Second, 
instead of all possible closed graphs, a regular lattice with 
edge-lengths rts 5b and 6 C was used. The dynamical so- 
lution inside the homogeneous region (that is, inside the 
horizon, where space is homogeneous but not static) was 
then obtained. An analytic continuation to the region 
outside the horizon shows that one can reduce the two 
free parameters by identifying the minimum area in the 
solution with the minimum area of LQG. The remaining 
unknown constant of the model, 5b, is the dimensionlcss 
polymeric parameter. This determines, with A m i n , the 
strength of deviations from the classical theory and must 
be constrained by experiment. 

The procedure to obtain the metric is in short the fol- 
lowing. 

i We define the Hamiltonian constraint replacing the 
homogeneous connection with the holonomies along 



ds 2 = -G{r)dt 2 + — — — - + H(r)dn 2 , 
F(r) 



G(r) 



F(r) = 



(r — r + )(r — r_)(r + r*) 2 



(r — r + ){r — r_)? 



H(r) = r 



(r + r,) 2 (r 4 + a 2) ' 

2 , a l 



(2) 



Here, r + = 2m and r_ = 



with dn 2 = dd 2 + sin 2 

2mP 2 are the two horizons, and r* = y/r+rZ = 2mP. P 

is the polymeric function P = (^/l + 7 2 — 1) /(vl + e 2 + 
1), with e< 1 being the product of the Immirzi parame- 
ter (7) and the polymeric parameter (5). With this, it is 
also P<1, such that r_ and are very close to r = 0. 
The area a is equal to A m ; n /87r, A m - m being the mini- 
mum area gap of LQG. Note that in the above metric, r 
is only asymptotically the usual radial coordinate since 
H(r) is not just r 2 . This choice of coordinates however 
has the advantage of easily revealing the properties of 
this metric as we will see. The ADM mass is the mass 
inferred by an observer at fiat asymptotic infinity; it is 
determined solely by the metric at asymptotic infinity. 
The parameter m in the solution is related to the mass 
M by M = m(l + P) 2 . 

If one now makes the coordinate transformation R — 
a /r with the rescaling t = tr 2 /a Q , and simultaneously 
substitutes R± = a /r^z, R* = a a /r* one finds that the 
metric in the new coordinates has the same form as in the 
old coordinates and thus exhibits a very compelling type 
of self-duality with dual radius r = ^JaZ. Looking at the 
angular part of the metric, one sees that this dual radius 
corresponds to a minimal possible surface element. It is 
then also clear that in the limit r — > 0, corresponding 
to R — > 00, the solution does not have a singularity, but 
instead has another asymptotically flat Schwarzschild re- 
gion. The dual mass observed from the dual observer in 
r w is m d = a a (l + P) 2 /AmP 2 . 



3 




FIG. 1: Penrose diagram of the regular static black hole so- 
lution with two asymptotically flat regions. Both horizons, 
located at r+ and r_, are marked in blue and red respec- 
tively. 



The causal diagram for this metric, shown in Fig [TJ 
then has two horizons and two pairs of asymptotically flat 
regions, A, A 1 and B,B', as opposed to one such pair in 
the standard case. In the region enclosed by the horizons, 
space- and timelikeness is interchanged. The horizon at 
r+ is a future horizon for observers in the asymptotically 
flat A, A' region and a past horizon for observers inside 
the two horizons. Similarly, the r_ horizon is a future 
horizon for observers inside the two horizons but a past 
horizon for observes in B, B' . If one computes the time 
it takes for a particle to reach r = 0, one finds that it 
takes infinitely long [f3|. The diagram shown in FigQ] 
can be analytically continued on the dotted horizons at 
the bottom and top. 

The metric in Eq. ([2]) is a solution of a quantum gravi- 
tationally corrected set of equations that, in the e, a — > 
limit, reproduce Einstein's field equations. However, due 
to these quantum corrections, the above metric is no 
longer a vacuum-solution to Einstein's field equations. 
Instead, if one computes the Einstein-tensor and sets it 
equal to a source term = SirT^, one obtains an 
effective quantum gravitational stress-energy-tensor T^ v . 
The exact expressions for the components of T are some- 
what unsightly and can be found in the appendix of [3f| . 
For our purposes it is here sufficient to note that the 
entries are not positive definite and violate the positive 
energy condition which is one of the assumptions for the 
singularity theorems. 

We can write the metric @ introducing the physical 



radius as radial coordinate. This is possible introducing 
the new radial coordinate x = \fH which varies in the 
range x e [y/2a , +oo [ 



x 



-f , x G [\/2a , +oo [. 



The metric assumes the following form 

2 



ds 2 



-G(r(x))dt 2 



1 



F(r(x)) \dx 



t\ da? 



(3) 



x 2 dn 2 , (4) 



where G, F are implicit functions of x. It is usual to de- 
fine g xx to be related to the mass inside a sphere of radius 
x and the gu component of the metric to be proportional 
to a dirty factor e 2 *( x ), since g t t ^ —l/g X x, 



g X x 
gu = 



2m(x) 



2m(x) 



=2*(x) 



and the dirty factor is 



2* 



(x) _ {r(x)+r*)\r\x)+al) 



(r*(s)-a2) a 

In the radial x coordinate the self-dual metric reads 



(5) 



(6) 



ds 



1 



2m(x) 



dt 2 



dx 2 



1 



2m(r) 



x 2 dQ 2 



The metric in this form satisfies the following Einstein's 
equations 



dm 
dx 



4irpx 



\ dG 2 (m(x) + 4ttP x x 3 ) 
G~dx~ ~ 
dP x = 
dx 



x(x — 2m(x)) 



^(P+Px) + ~(P±-Px), (7) 

Lr ax x 



where energy density and pressures are defined by T M „ = 
G M i//87r = Diag(— p, P x , P±, P±) using the coordinates 
(t,x,9,4>) and the metric (@|. The first equation in ([7]) 
shows that m(x) is the mass inside a shell of radius x and 
from the first equation in ([5]) we can extract the function 
m(x) (r = r(x)), 



m(x) 



1 a 2 



+ r 



(r 4 — a 2 ) 2 (r — r_)(r — r + ) 
(a2+r 4 ) 2 (r + n) 2 



where r = r(x) is defined implicitly by the relation Q. 
The mass m(x) above tends to the ADM mass for r — > 
+oo and to the dual ADM mass for r — > (or x —> +oo), 



m(x) 



m(l + P) 2 for x — > +oo , r > ^/a^, 



4mP 2 



for x 



-oo , r < Joo. 



(8) 
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The exact relation between the physical radial coordinate 
x and r is 



T 2 \l x^— 4 a? r ( 

— - ~ J for r < V^°> 



(9) 



for r > xfaZ, 



III. VAIDYA COLLAPSE 

We will proceed by combining the static metric with a 
radially ingoing null-dust, such that we obtain a dynam- 
ical space-time for a black hole formed from such dust. 
In the present model this process, usually described by 
the Vaidya metric jifj, will have corrections, negligible 
in the asymptotic region, but crucial to avoid the for- 
mation of a singularity in the strong-curvature region. 
The metric constructed this way in the following is not 
a strict solution of the LQC minusuperspace equations 
in the vacuum. In other words, keeping in mind what 
we said in the previous section, it is a solution of the 
Einstein equations with an effective energy tensor which 
depends also on the time variation of the mass. 

We use the radial coordinate r of the previous section 
and we start by making a coordinate transformation and 
rewrite the static space-time in terms of the ingoing null- 
coordinate v. It is defined by the relation dv = dt + 
dr/y/F{r)G{r), which can be solved to obtain an explicit 
expression for v. The metric then takes the form 



ds 2 



-G{r)dv 2 



I G(r) 
F(r) 



drdv + H(r)dfL 2 . (10) 



Now we allow the mass m in the static solution to depend 
on the advanced time, m — > m(v). Thereby, we will 
assume the mass is zero before an initial value v a and 
that the mass stops increasing at v^. We can then, as 
before, use the Einstein equations G = 8nT to obtain 
the effective quantum gravitational stress-energy tensor 
T. T v v and T r r do not change when m(v) is no longer 
constant. The transverse pressure T e g = Tt however has 
an additional term 



T e e (m(v))=T s e (m) 



Pr 2 m'(v) 



27r(r + 2m(v)P) 4 



(11) 



where m' = dm/dv. Because of the ingoing radiation, the 
stress-energy-tensor now also has an additional non-zero 
component, T r v , which describes radially ingoing energy 
flux 



This metric reduces to the Vaidya solutions at large ra- 
dius, or for e — > 0, a — > 0. However, in the usual Vaidya 
solutions, the ingoing radiation creates a central singular- 
ity. But as we see here, with the quantum gravitational 
correction, the center remains regular. 

We note that the ingoing energy flux has two zeros, one 
at r = r*(v) and one at r — y/aZ, and is negative between 
these. What happens is that the quantum gravitational 
correction works against the ingoing flux by making a 
negative contribution until the effective flux has dropped 
to zero at whatever is larger, the horizon's geometric 
mean r* or the location of the dual radius r = 
The flux then remains dominated by the quantum grav- 
itational effects, avoiding a collapse, until it has passed 
r* and the dual radius where it quickly approaches what 
looks like an outgoing energy flux to the observer in the 
second asymptotic region. 



IV. A SCALAR FIELD ON THE LBH 
BACKGROUND 



The wave-equation for a massless scalar field in a gen- 
eral spherically symmetric curved space-time reads 



1 



=0 M (g 



gd v 9) = o, 



(13) 



where <f> = $>(r,8,<j),i). Inserting the metric of the self- 
dual black hole we obtain the following differential equa- 
tion 



H(r) 2 



.<9 2 $ 



G{r)F'(r)& 

+ cote— + 



- 2G <'-> («r 

-F(r) \H{r)G'{r)<S>' + 2G(r) (h'(t)<P' + H{r)<$> 



(14) 



= 0, 



where a dash indicates a partial derivative with respect 
to r. Making use of spherical symmetry and time- 
translation invariance, we write the scalar field as 



$M,<M) :=T(t)<p(r)Y(0,, 



(15) 



omitting the indexes I, m in the spherical harmonic func- 
tions Yi m (0,(j>). The standard method of separation of 
variables allows us to split Eq. p3| in three equations, 
one depending on the r coordinate, one on the t coor- 
dinate and the remaining one depending on the angular 
variables 0, <h. 



2(1 + P) 2 r 4 (r 4 - a 2 )(r - r*(v))m'(v) 
(a 2 + r 4 ) 2 (r + r t (v)) 3 



(12) 



Notice that also in the dynamical case, trappin g horizons 
still occur where g rr — F(r, v) vanishes [42,|43|], so we can 
continue to use the notation from the static case just that 
r±(v) and r»(u) are now functions of v. 



GF 9 f HVGF ^ 
dr 



H dr 

do 2 

dt 2 



d 

cote— + csc 2 e' 
de 

T(t) = -cu 2 T(t), 



r,K l + 1) 

Y{6,4>) 



-K 2 Y{e,cj>), 



(16) 
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where K 2 = 1(1 + 1). To further simplify this expres- 
sion we rewrite it by use of the tortoise coordinate r* 
implicitly defined by 



dr* 
dr 



1 



IGF 



(17) 



Integration yields the new radial tortoise coordinate 

(18) 
log \r -r+\ . 



2 (r_ + r + ) 

Q Q „ 2 „2 l0 S( r ) 



.(r+ - r_) 



log |r 



_(r+ - r_) 



Further introducing the new radial field <p(r) 
t/}(r)/y/H the radial equation ([T6|) simplifies to 



(9 



9r 



F(r) = 



- +LU 2 -V(r(r*)) 



V>(r) = 0, 



(19) 



1 / GF 

IT 




Inserting the metric of the self-dual black hole we finally 
obtain 



V(r) = 
2 



(r — r_)(r — r + ) 
(r 4 + a 2 ^ 



r ( (K 2 - 2) r + r_ + r+) + 2K 2 rr* + K 2 r 2 ) 
+2a 2 r 4 ((K 2 + 5) r 2 + 2K 2 rr* + K 2 r 2 - 5r(r_ + r+) 
+5r_r + ) + r 8 (A" 2 (r + r») 2 + r(r_ + r+) - 2r_r+) 



The potential V(r) is zero at r = r + and r_ as for the 
classical Reissner-Nordstrom black hole. We therefore 
can follow the same analysis as for this case, approxi- 
mating V(r(r*)) near the horizons via 

V(r*) oc e 2K+r , for r — > r + or r* — ^ — oo , 
V(r*) cx e~ 2K ~ r , for r — > r_ or r* — » +oo, 
T^(r*) — S- , for r -> or r — > +oo. (20) 

A crucial quantity to study black hole evaporation is the 
surface gravity 



F(G'(r)f 
AG 



(21) 



where ' denotes the derivative respect to the radial coor- 
dinare r and = (1,0,0,0) is a timelike Killing vector 
in r > r+ and r < r_ but space-like in r_ < ?* < 7'+. For 
the metric © we find the following values 



4m 3 P 4 (l - P 2 ) 



4m 3 (l - P 2 ) 



, (22) 



16to 4 + a 2 

for the surface gravity on the inner and outer horizons. 



V. PARTICLE CREATION 

We now work out the particle production in the back- 
ground geometry describing the formation of a LBH. The 
spacetime associated to the gravitational collapse to form 
a black hole is not everywhere stationary and then we ex- 
pect particle creation. The particle creation phenomenon 
is due to the non stationary gravitational collapse. How- 
ever the spacetime will be stationary at late time and 
then particle creation is just a transient phenomenon and 
we can forget all the details related to the gravitational 
collapse. We introduce the Hawking effect in the sim- 
plest possible scenario a la Vaidya explained in section 
piip . The basic tools to evaluate the particle production 
are the Bogulibov transformations needed to connect the 
positive frequency modes of the field between the initial 
and final stationary regions . The field can in fact be 
expanded in the initial stationary region as 



a^ fu + <\j T f u 



(23) 



(the use of discrete modes from the continous one can be 
implemented smearing with suitable wave packets 0, [20l 
l45j ) or in the final one as 



(24) 



where a™,a° ut are the ladder operators verifying the 
usual commutation relations and f w and p^ are the solu- 
tions of (|TT))) in the initial and final regions respectively. 

We define the in Fock space with the natural time v at 
I~ . The positive frequency modes which are solution of 
(p~9|) at infinity (I - ) are 



fu(r, v) 



-iuj(r* -\-t) 



(25) 



where v = t + r* and they obey the scalar product 

(/o))/w') = — (/ul/u 1 ) = 

-ij do Hd£l(f u d v fo 6(w-u/) (26) 

and (fuift)') ~ 0- We define also the out Fock space at 
I + associated with the natural time parameter u. The 
positive frequency modes are 



Pu(r> u ) 



)-* -ft) 



47Tv / W 1 /ff AlTy/cHy/H' 



(27) 



where u = t — r* and the outgoing modes obey the nor- 
malization condition 

(p u ,Pu') = -(pZ,pI') = 

-i / duHdQipud-aPu* - pl,d u p u ) = 5(u> - u') (28) 
Ji+ 
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and (jPu-iPu>) — 0- The modes f u and p u are solutions of 
(|19[) at I~ and / + respectively; we can also approximate 
H w r 2 in these the asymptotic regions. We decided 
about / + as Cauchy surface but it is not properly correct. 
We must include the future event horizon H + to have a 
complete Cauchy surface: I + 1J H + . The modes p u are 
not complete and we have to add those that cross the 
future horizon H + . However we do not need them to 
evaluate the particle production at I + because the result 
is insensitive to the ingoing modes. The next step is to 
calculate the Bogoliubov coefficients relating the ingoing 
and outgoing basis solutions f u and p u . Assuming the 
frequency to be discrete 

Plo ^ ^ A UUJ ' f u ' -\- B ULlJ 'f u i^ (^9) 

and 

A uu > = (p u ,fu>) & B uu ' = —(Pu,fu>) (30) 

which satisfy the following matrix relations again assum- 
ing discrete values for the frequency, 

AA f - BB ] = 1, 

AB T - BA T = (31) 
and the matrix elements are 

A^ = -i^ dvHdn(p u d v f* - f*,d vPu ), (32) 

where, for mathematical convenience we have chosen I~ 
as the Cauchy surface to calculate the scalar product that 
is insensitive to this choice. The Bogulibov coefficients A 
and B can also be used to expand one of the two sets of 
creation and annihilation operators in terms of the other, 

^ = J2 A ^<' t+B ^ a °^ ( 33 ) 



< ut = J2 A ^' a -'- B ^'<' ( 34 ) 

If any of the B^^ 1 are non zero the particle content of 
the vacuum state at I~ (which we indicate with the ket 
|in)) respect to the Fock space at I + is non trivial, 

(in|A^ + |in)=£|5_f • (35) 

where is the particle number operator at frequency 
u at J + . In contrast, if all the coefficient B^,^ are equal 
to zero, the first of the relations ([31]) reduces to AA^ = 1 
and then the positive frequency mode basis f u and p^ are 
related by a unitary transformation and the annihilation 
operators (|33|) and (|34|) define the same vacuum. 

To evaluate the products (|3"D)) we need to know the 
behavior of the modes p u at I~ . For this propose we 




FIG. 2: Penrose diagram for the Vaidya collapse and slowly 
evaporation. The Penrose diagram for the whole collapse and 
evaporation process will be given in the next figure. 

consider a geometric optic approximation in which the 
massless particle world- line is a null ray, 7, of constant 
u and we trace this ray backwards in time from I + until 
I~ (see Fig 12]). The later it reaches / + , the closer it must 
approach H + . The ray 7 is one of the rays whose limit as 
t — > +00 is a null generator jh of H + . To specify 7 we 
can then use its affine distance from along an ingoing 
null geodesic through H + . This can be easly found using 
the Kruscal-Type coordinates in the re gion outside the 
horizon r + = 2m. These are defined by [ll( 

U+ = -— e- K +" , V + = —e K + v , (36) 

where k + is the surface gravity calculated in (f2"2"]l. The 
affine parameter in the ingoing null geodesic closer to H + 
is U + = — e and so using the first of (|36)) we find 

it = log e + const. (37) 

K+ 

substituting the previous expression in the solution (|27p , 
we can see that it oscillates rapidly at later times t and 
this justifies the geometric optics approximation. We 
need to match p^, with the solution of the Klein-Gordon 
equation near I~ . In the geometric optic approximation 
we just parallel-transport the vector n a , tangent to the 
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null geodetic which is ingoing at H + , and l a which is the 
null generator of H + , back to l~ along the continuation 
of 7jj . We call vq the point where the continuation meets 
/ _ then the continuation of the ray 7, along the outgoing 
null geodesic, meets I~ at v = Vq — e so 



■ \og(vo — v) 



for V < V , 



Pu = for v > v , (38) 

where vq is the latter time at which the field can reach 
infinity without entering in the black hole. Next step is 
the calculation of the Bogolubov matrix A. Introducing 
(l38l) and (ESI in (El we find 



A , - 
B. 



/ K + ( UJ 
2ir\/ojui' V K + 



0' ^-^-OJ,— UJ' 1 



(39) 



where we define Vq = 0. The Bogoliubov coefficient Aj.u' 
is the fourier transform of a function that vanishes for 
v > vq than it is analytic in the lower half of the complex 
oj' plane. It has a logarithmic branch point in a/ = 
then the branch cut extends into the upper half plane. 
Therefore, we have the following relation between the 
coefficients A and B 

|4^| = e^|fl„ lU '|. (40) 
Finally introducing (|40| in the first relation (|31[) we find 

= e ^ +u 'V K + - ll (Bflt) UiU)i (4i) 

Taking uj = us' the number of particles (|3"5j) in the w th 
mode at 1+ is 



(in|iVi + |in) = ( J BBt )[ , 



(42) 



The result (|42|) coincides with the Planck distribution of 
thermal radiation for bosons at the temperature Tbh = 

K+/27T. 

Evaporation time. The evaporation proceeds through 
the Hawking emission at r+, and the black hole's 
Bckcnstcin-Hawking temperature, given in terms of the 
surface gravity k by T BH = k/2it, yields [l3j 



T B H(m) 



(2to) 3 (1-P 2 ) 
47r[(2m) 4 + a 2 ] 



(43) 



This temperature coincides with the Hawking temper- 
ature in the limit of large masses but goes to zero for 
??i -> (FigE}. 

The luminosity can be estimated by use of the Stefan- 
Boltzmann law L(m) = aA B (m)Tg H (m), where (for a 
single masslcss field with two degrees of freedom) a = 




FIG. 3: Plot of the temperature versus the loop black hole 
mass. The dashed line represents the classical temperature 
T = 1/8™. 



tt 2 /60, and A H (m) = 47r[(2m) 2 + a 2 /(2m) 2 ] is the surface 
area of the horizon. Inserting the temperature, we obtain 



L(m) 



I6m 10 a(l-P 2 ) 4 
7r 3 (a 2 + 16m 4 ) 3 



(44) 



The mass loss of the black hole is given by the following 
equation for the black hole mass 



dm 
dv 



-L{m) 



(45) 



and we can integrate its inverse to obtain the mass func- 
tion m(v). The result of this integration with initial con- 
dition m(v = vq) = toq is 



Az; 



{5a 6 a + 432o^to 4 + 34560a 2 m 8 - 61440to 12 )tt 3 



(5a 6 + 432a>o= 



720to 9 (1 - P 2 ) 4 a 
- 34560a 2 TO i - 61440toJ 2 )tt 3 



720to '(1 - P 2 ) 4 a 



(46) 



In the limit to — > this expression becomes Av w 
aQ7r 3 /(144TO 9 (l — P 2 ) 4 a), and one thus concludes that 
the black hole needs an infinite amount of time to com- 
pletely evaporate. In the complete evaporation process, 
we neglected the backreaction for any value of the mass 
because at v ~ +00, when to < top, dm/dv << to and 
then, contrary to the classical case, such approximation 
is valid also in the final stages of evaporation. 



VI. COLLAPSE AND EVAPORATION 

We are now ready to combine the black hole formation 
and evaporation. As in section UTTT we divide space-time 
into regions of advanced time. We start with empty space 
before v a then the mass increase from v a to Vb since the 
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gravitational collapse. For astrophysical black holes this 
evaporation will proceed very slowly, and m remains con- 
stant to good accuracy at Too, but at some later time, v c , 
Hawking radiation becomes relevant and to decreases un- 
til it reaches zero again in an infinity time as we have seen 
in the previous section. 

We thus have the partition — oo < v a < v b < v c < oo 
with 



Vv € (-00, v a ) 
\fv e {v a ,v b ) 

Vu e (v b , v c ) 
Vw e (v Cl +00) 
for v — > +00 



m(v) = 0, 
m'(v) > 0, 
m(v) = mo, 
m'(v) < 0, 
m(v) — > 0. 



(47) 



The mass would immediately start to drop without in- 
coming energy flux and thus v a — Vb, but stretching this 
region out will be more illuminating to clearly depict the 
long time during which the hole is quasistable. 

To describe the Hawking-radiation we will consider the 
creation of (massless) particles on the horizon such that 
locally energy is conserved. We then have an ingoing 
radiation with negative energy balanced by outgoing ra- 
diation of positive energy. Both fluxes originate at the 
horizon and have the same mass profile which is given by 
the Hawking temperature. The area with ingoing neg- 
ative density is again described by an ingoing Vaidya 
solution, while the one with outgoing positive density is 
described by an outgoing Vaidya solution. 

The outgoing Vaidya solution has a mass-profile that 
depends on the retarded time u instead of v and the 
mass decreases instead of i ncreases. T he retarded time is 
defined by du = dt — drj y 'F(r)G(r). After a coordinate 
transformation, the metric reads 



ds 2 



-G(r, u)du 2 



F{r,u) 



dudr + H{r)dfl 2 , (48) 



where F(r,u) and G(r, u) have the same form as in the 
static case but with m replaced by a function m(u). 
We fix the zero point of the retarded time u so that 
r = r + corresponds to u c = v c . Then there is a static 
region with total mass mo for v > v Cl u < u c . Note that 
since the spacetime described here has neither a singu- 
larity nor an event horizon, we can consider pair creation 
to happen directly at the trapping horizon instead of at 
a different timelike hypcrsurface outside the horizon, as 
done in 44|. We have in this way further partitioned 
spacetime in regions, broken down by retarded time: 



Vu < u c : m(u) — too , 
Vu > u c : m'{u) < . 



(49) 



Now that we have all parts together, let us explain the 
complete dynamics as depicted in the resulting causal 
diagram FigJU 

In the region v < v a we have a flat and empty region, 
described by a piece of Minkowski-space. For all times 




FIG. 4: Penrose diagram for the formation and evaporation 
of the regular black hole metric. The red and dark blue solid 
lines depict the two trapping horizons r_ and r+. The brown, 
dotted line is the curve of r = sJaZ and the brown, long dashed 
one is r*. The light blue arrows represent positive energy flux, 
the magenta arrows negative energy flux. 



v > v a , the inner and outer trapping horizons are present. 
These horizons join smoothly at r — in an infinite time 
and enclose a non-compact region of trapped surfaces. 

A black hole begins to form at v = V a from null dust 
which has collapsed completely at v = v b to a static state 
with mass mo. It begins to evaporate at v = v c , and the 
complete evaporation takes an infinite amount of time. 
The observer at / + sees particle emission set in at some 
retarded time u c . The region with v > v c is then divided 
into a static region for u < u c , and the dynamic Vaidya 
region for u > u c , which is further subdivided into an 
ingoing and an outgoing part. 

As previously mentioned, the radially ingoing flux 
(light blue arrows) in the collapse region is not positive 
everywhere due to the quantum gravitational contribu- 
tion. It has a flipped sign in the area between r* (black 
short dashed curve) and r = ^fdZ (brown dotted curve) 
which is grey shaded in the figure. Likewise, the ingoing 
negative flux during evaporation (magenta arrows) has 
another such region with flipped sign. It is in this region, 
between the two horizon's geometric mean value r% and 
the dual radius corresponding to the minimal area, that 
the quantum gravitational corrections noticeably modify 
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the classical and semi-classical case, first by preventing 
the formation of a singularity, and then by decreasing the 
black hole's temperature towards zero in an infinity time. 

VII. UNITARITY PROBLEM 

In this section following Wald book [45[ , we resume the 
main features of the loss of information and unitarity of 
classical black holes. The evolution of a physical system, 
represented by a quantum state |in) in an initial Cauchy 
surface in Minkowski space, is given by a unitary op- 
erator that maps it into a final state |/) in a final Cauchy 
surface Let's see what happens when a classical black 
hole is present: The Cauchy surface £/ can split into two 
surfaces: £/ = Ej nt 1J E ut of which the first lies inside 
the black hole and the second lies outside. We can also 
think as limiting case Sj nt = 11^+ and E ut = £ /+ . Now 
if |in) = |0)/- we know from [45j that we can formally 
write 

U\Q)j- = 1] I £ e-^\Ni)- mt g> |iVi) out 1 , (50) 
% Ijv=o J 

with U given by a unitary transformation relating the 
vacuum of the initial Fock space to the N particle states 
of the final Fock space. The form of the right hand side 
is due to the presence of a bifurcated Killing horizon and 
the consequent different notions of time on the initial and 
final Cauchy surfaces. Note that the expression ([50]) is 
formal because the state on the right hand side fails to 
be normalizablc if the condition 

Tr(B*B) < oo (51) 

is not fulfilled j4j|. This condition implies that the total 
number of particles produced in the transition between 
the initial and final state is finite. If it is not satisfield 
the initial and final Fock spaces are unitary inequivalent, 
neverthless the expression ([5U|) can still be used to de- 
scribe approximated states (see HBj]). Note in fact that 
in both the usual Unruh and Hawking effects, this is the 
case. 

Now the presence of the horizon obliges the exter- 
nal observer to trace on the internal degrees of freedom, 
transforming the initial pure state |in) in a density ma- 
trix p. This is the information loss. Thus it is the horizon 
that causes the presence of a mixed state with an asso- 
ciated temperature, as a consequence of our ignorance of 
the complete system. The state that the observer at I + 
sees will not only be mixed but also made of uncorre- 
lated radiation (stochastic thermal radiation). The im- 
portant point however is that the external observer loses 
the correlation between the two regions due to the trac- 
ing operation: but the correlations do exist!. Altought 
the information is lost we can claim that the quantum 
unitary evolution is preserved at a fundamental level. If 
we were able to look inside the horizon we would see the 
correlations and the state would continue to be pure. The 



i+ 



r = 




i 



FIG. 5: Hawking scenario 



particles propagating to infinity are strongly correlated 
with particles that enter the black hole at early times. 
The presence of a density matrix to describe the state of 
the field in the exterior region from the algebraic point 
of view simply corresponds to the fact that the restricion 
of the state of the field to the subalgebra associated to 
the int D(S out ) (where int means interior of a set and D 
is the domain of dependence, see |45[) leads to a mixed 
state and this is due to the fact that the '"domain of de- 
terminacy"' of int D(Y, ont ) is not the entire spacctime. 
The breakdown of the unitarity arises from the presence 
of a true singularity in the case of complete evaporation 
see Fig(5j 

The disappearence of the black hole would completly 
remove from the spacetime the correlations (that could 
have preserved purity) hidden by the Horizon. As a con- 
sequence, on a surface £^ in the final spacetime, after the 
evaporation, we would have really a mixed state also at 
a fundamental level. The conclusion is that not only we 
have information loss in the process |in) — > p, but also 
that this process is not Unitary. The problem is that 
now at late time, the entire algebraic state of the field 
is mixed. In algebraic terms we can look at the whole 
process as an evolution of the state in a time labeled by 
the Cauchy surfaces. In the figure ([5]) we are looking 
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at the evolution from the surface S/ to the surface T/j. 
Now the domain of determinacy of the hypersurface S f 
is the entire spacetime (see |45(). On the other hand 
int Z?(Sj) includes the entire future of the spacetime 
(after the evaporation), but its domain of determinacy 
is not the entire spacetime, since it does not include the 
black hole region. The evolution corresponds then to the 
restriction of a pure state u> from the Weyl algebra of the 
entire spacetime to the subalgebra associated to the re- 
gion int -D(E'j) i.e evolution from a pure to a mixed state. 
The entire problem can be summarized saying that we are 
evolving from a Cauchy surface to a surface that fails to 
be a Cauchy surface for the whole spacetime. In the next 
section we will see in which sense these problems can be 
cured in our model. 



VIII. UNITARITY RESTORED 

We have shown in this paper that LBHs evaporate 
emitting particle in thermodynamical equilibrium at the 
temperature Tbh = k + /2tt. The black hole evaporation 
time is infinite (|46[) . As we have seen in the previous 
section, the black hole paradox presents essentialy two 
features: the loss of information, strictly linked to the 
existence of an horizon, and the not unitary evolution, 
that is essentially related to the presence of a singular- 
ity. The effective spacetime that we have described in 
the previous sections seems to be a good candidate to 
solve both problems. The spacetime is in fact singularity 
free and the Bogoulibov coefficients evolve with the mass 
allowing a unitary tranformation between the initial and 
final vacuum. Let's see in detail the behaviour of these 
coefficients: for small but non zero black hole mass we 
have thermal radiation but when the mass goes to zero, 
in an infinite amount of time, the unitarity is restored. 
More concretely looking at the surface gravity k + , we see 
the mass goes to zero (in an infinite time) together with 
the temperature. This is the crucial difference with the 
classical explosive case where the temperature goes to in- 
finity when the mass reduces to zero. From the Planck 
spectrum in (|42p we deduce that in an infinite amount of 
time 



lim {BB% iUI = Vw, 

m— >0 



(52) 



because n — > for m — > 0. Since {BB^) UtU1 is positive 
semi-define, this vanish iff i? = 0. Let us now recall 
another useful relation between the initial and final Fock 
spaces. Using the Bogoliubov transformation between 
creation and annihilation operators respectively at I~ 
and J + , we can find the relation between the vacuum 
state | in) (at I") and the vacuum state |out) (at I + ), 



| in) = (out | in) 
where the matrix V is 



| out), 



-B* ,„A; 



(53) 



(54) 



Now for m 
in itself 



0, B — >• and then the initial state evolves 



| out). 



(55) 



We conclude the unitarity is restored in an infinite 
amount of time. Since B = the positive frequency 
mode basis f u and are related by the unitary transfor- 
mation A as immediate consequence of the first relation 
in (|3"Tj) . Note also that due to this behaviour the |out) 
state, unlike in the usual treatment, can be related to the 
initial vacuum state by an exact unitary transformation 
because the necessary condition (|5l"j) is satisfied. 

To test the radiation emitted is thermal, for non zero 
black hole mass, we can calculate the probabilities of 
emitting different numbers of particles. For instance, we 
can calculate 120] 



(m\Ni + Nt + |in) = 1 V -. NO >-. 



(56) 



which agrees with a thermal distribution. In a similar 
way we can find all the higher moments that coincide 
with the thermal probability 



P(N U ) = (1 _ e -W«+) e -2TJVo;/K + _ 



(57) 



to emit N particles in the mode u>. If we wait for an 
infinity amount of time the mass goes to zero together 
with the thermal emission probability 

lim P(N U ^ 0) = , lim P(N„ = 0) = 1 (58) 

and the particle emission stops after an infinite amount of 
time. Using ([57]) the expectation value of the out particle 
number operator reads 



(in|7Vi + |in) 



+ 00 

E 

N=0 



NP(N U ). 



(59) 



The emission probability yet meets the normalization to 
one in the zero mass limit. 

To check the thermal properties of the radiation pro- 
duced by the black hole we calculate 

(in|iVi + ^|in) = (in|7Vi + |in) (in|^|in), 

this result shows the absence of correlations between dif- 
ferent modes as typical of the thermal radiation, then the 
quantum state at / + is exactly described by a thermal 
density matrix 



N=0 



\N U )(N U \, (60) 



where |iV w ) is the state at I + with TV particles in the 
mode iij. We see that any measure at is described 
by the density matrix p and the von Newmann entropy 
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SVad = — Trplogp associated to the frequency mode w 
reads 



QUI 



rad 



2lTUJ I K + 

2ttui/k + Y 



log 1 



-27TLJ / K,+ 



(61) 



The annihilation operators for the modes entering the 
black hole region are then 



E*. ,t€L] , 



mj 



(68) 



and the total entropy is SVad = J °° dwS^ = 7tk + /6. 
When the balck hole mass goes to zero the density matrix 
reduces to 



hm /9 = ni0w)(0w| = |in)(in|, (62) 
and the state is pure. The process in its totality is 



p{m > 0) -> |in)(in|, 



(63) 



showing the initial pure state at J~ evolves to itself at 
/ + . The radiation at I + is no more thermal in the limit 
m — > and the hnal state is not a mixed state; the 
complete evolution is unitary. The entropy (|5T|) goes to 
zero. 

During its evaporation the black hole is in thermal 
equilibrium with the radiation outside for each finite 
value of the mass. When the black hole mass reduces 
to zero (in an infinity amount of time) the matter out- 
side is no more in a thermal state but in a pure state. 
Thermal radiation is present outside the black hole only 
for finite values of the black hole mass but when m w 
the mixed state collapses to a pure state. 



Because we have ignored this sector we have lost all the 
possible correlations with the quanta entering into the 
black hole. If we properly take into account those modes, 
the purity of the |in) state is restored. The state ([53")) is 
now 



= (H 



me 2 



|fin), 



(69) 



with the index i, j variyng over all the modes on the final 
Cauchy surface (int, and out), and |fin) = |0)j„t ® |0) OU f. 
Following [20] we find an explicit relation between |in) 
and the TV-particles states with frequency uj at / + and 
H + , respectively |A^ + ), \N, 



n 



E 

N=0 



1-e - + > e - + \Nl )®\N? ). (70) 



This relation shows we have an independent emission 
in frequencies of quantum entangled states representing 
outgoing and ingoing radiation. Note that for m 7^ 
the thermal density matrix is recovered evaluating the 
expectation values of any operator O at / + , for example 



A. Role of the observer and correlations 

In the previous section we have restriced our attention 
to the surface I + , now we condider also the previously 
neglected modes. When the black hole mass is non zero 
the pre state |in) could appear mixed but it is still a pure 
state since we have to add the modes crossing the event 
horizon H + . A proper expansion of the field is then 



E 



a out V 



int i 



n u ' T ir 



(64) 



where aj^' are the ladder operators for the incoming par- 
ticles at the future horizon and h u are the corresponding 
modes. The modes h u are not uniquely defined because 
there is no natural time parameter on H + . A simple 
expression for these modes [2(J EH can be obtained re- 
versing the sign of vq — v and us in (|3"5)l . With this choice 
the additional Bogulibov coefficients G, E defined by 



(65) 



where 



Gww' = (/iw,/w') & -Ew = ~{h u , ft,') (66) 
have simple relations [2(j with the coefficients A and B, 



G* 



(67) 



<in|JVi + |in) =Tr(pA /+ ) 



(71) 



Tr 



E 

N=0 



e~^\N^){Ni + \N I+ 



In Fig|4] we can see that the evolution is unitary at 
a fundamental level, in fact the absence of singularity 
implies that any Cauchy surface of the spacctime is a 
Cauchy surface for the whole spacctime and the Hawking 
scenario of Figj5]is avoided. If we use the Cauchy surfaces 
E as a "time label" for the spacetime describing the whole 
evaporation, we don't find any discontinuity in contrast 
to the Hawking scenario. The presence of a continuum 
shrinking of the mass, allows us to describe the state 
(|50|) or equivalcntly (|70[) in terms of a mass dependent 
splitting between internal and external spacctime of the 
final Cauchy surface Ef = Ej™. [j E™ t : 



moh- = u E 



■\N„)Z t ®\N u )l 



(72) 



,N=0 



which is based on the decoposition of the final one par- 
ticle Hilbert space in terms of the direct sum of two one 
particle Hilbert spaces, one internal and the other ex- 
ternal to the horizon, = H™ t © ff^ t . When m 
is zero the Hilbert space of solutions reduces to H^ n = 
H I=0 8 H 1+ . For m -> the relation (J72J) tends to 



^|0) 7 - -> |0) int ® |0) 

c 



(73) 
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FIG. 6: Role of the observer and correlations 

but also B — > and as a consequence of (|67l ) E — > 0. In 
this limit then the initial vacuum reduces to 

|0) in - |0)out ® |0) int - |0)fin (74) 

and the density matrix to the pure state (|G2 j) . 

We can better understand the properties of the radi- 
ation introducing an abserver in the classical region of 
the dynamical spacetime describing the collapse. The 
situation is depicted in Figure The Observer 0-2 at 
u << u c will not see any remarkable radiation because 
he is too far from the event horizon. Any Observer O3 for 
— 00 > u > u c will see thermal radiation as an incoming 
flux of particles from the region where the matter col- 
lapsed. The essential point is the presence of an horizon: 
the pure state is perceived as a flux of particles because 
the observer is obliged to trace over the internal degrees 
of freedom, losing the information about the particles in- 
side the black hole. The Observers O3 will experience 
a thermal particle bath of increasing or decreasing tem- 
perature depending on the stage of the evaporation; the 
temperature is given by (|43|) . Finally for the observer 0\ 
(which follows a trajectory with constant radial coordi- 
nate) at i + the initial vacuum state evolves in itself as 
evident from ([72"]) when n — > 0. The 0\ observer sees a 
gas of particles of increasing or decreasing temperature 
moving in time towards i + . During the evaporation the 
state is not pure because the trace on the modes that 
cross the future horizon H + but when the mass goes to 



zero unitarity is restored. It makes sense to speak about 
density matrix for any non zero value of the black hole 
mass but for m « such operator reduces to a pure state 
in a natural way. We wish to mention the fact that the 
observer in i + is in causal contact with both the radia- 
tion outcoming towards / + and the one ingoing to r = 0, 
this fact in principle allows him to see the correlations 
between these modes and in this sense he could detect 
the previously hidden correlations needed to restore pu- 
rity at all the times (a detailed study of of these corre- 
lations is currenlty under investigation [47j) Note that 
any observer not located at / + would see a Temperature 
To = k/2-kx that differs from the temperature T by the 
ratio of "Killing time" to "proper time" for that observer, 
where x is the redshift factor given by the square of the 
norm of the Killing vector. 

We conclude this section with a comment. We can 
compare the result in this section with what happens 
in the Rindler spacetime. The Unruh temperature of the 
radiation is Tjj — g/2ir and it goes to zero when the accel- 
eration g vanishes. The Bogoliubov transformation also 
goes to zero with the acceleration. For the LBH when 
the mass is smaller the Planck mass the behavior is the 
same, the temperature goes to zero with the decreasing 
black hole mass. In other words the limits of zero ac- 
celeration and empty space coincide leaving a pure state 
and the full process is unitary. 



IX. CONCLUSIONS 

We have calculated the particle creation by a loop 
black hole that is entirely singularity-free. The analy- 
sis follows the classical one for the Schwarzschild space- 
time but with the crucial difference residing in the sur- 
face gravity of the black hole. The classical surface grav- 
ity (proportional to the temperature) diverges when the 
black hole mass goes to zero. For the LBH instead it goes 
to zero with the mass but in an infinite time. The approx- 
imation we used in this paper is that the matter fields (in 
particular we concentrated on the massless scalar field) 
obey the usual wave equations with the Minkowski met- 
ric replaced with a black hole spacetime metric g^ v which 
is solution of Einstein equations with an effective energy- 
tensor builds up a negative contribution that violates the 
positive energy condition and prevents the formation of 
a singularity. We went trough the details of the Hawking 
calculation for a massless scalar field that propagates in 
the given fixed black hole background and we found the 
Bogoliubov transformation as function of the new surface 
gravity. In an infinity amount of time the surface gravity 
goes to zero together with the Bogoliubov transformation 
and the whole process, collapse and complete evaporation, 
is unitary. The model we have studied in this paper is 
very simple and the analysis very similar to the classi- 
cal one but the quantum gravity properties of the metric 
imply that the vacuum pure state at I~ evolves in itself 
at I + . We designed also the causal diagram for the com- 
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plete process of collapse and evaporation. The value of 
the scenario studied here is that it provides a concrete, 
calculable, model for how quantum gravitational effects 
alter the black hole spacetime and the particle flux to I + . 
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